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Abstract 

A collider signal with a stable gravitino of O(10)eV mass at the International 
Linear Collider (ILC) experiment is investigated. Such a light gravitino is gen- 
erally predicted in the low-scale gauge mediation scenario of the supersymmetry 
breaking. We particularly focus on the case that the next lightest supersymmet- 
ric particle is stau, which eventually decays into a gravitino and a r-lepton. With 
such a small gravitino mass, the lifetime of the stau is 10~^^-10~^^sec, and the 
produced stau decays before reaching the first layer of the inner detector of the 
ILC. It is shown, however, that the lifetime can be determined from the distri- 
bution of the impact parameter, which is obtained by observing charged tracks 
caused by decay products of the r-lepton. This measurement also enables us to 
estimate the mass of the gravitino and determine the scale of the supersymme- 
try breaking. Based on a simulation study, we found that the lifetime can be 
measured when it is longer than ~ 10~^^sec and the stau mass is about lOOGeV. 



1 Introduction 



The Large Hadron Collider experiment (LHC) is now operating in order to clarify 
new physics beyond the standard model (SM). Though positive signals of the new 
physics are, unfortunately, not detected so far, those are expected to be found in 
near future, because the hierarchy problem of the SM strongly suggests that the new 
physics should be at the TeV scale or below. On the other hand, many new physics 
models have been theoretically proposed so far. Among those, the supersymmetric 
model [2] is the most attractive one, because it guarantees the stability of the Higgs 
mass against radiative corrections and gives a clue to solve the hierarchy problem of the 
SM. In addition, the supersymmetry (SUSY) plays a crucial role to realize the grand 
unification of known gauge interactions of the SM at a certain high energy scale. 

Details of the supersymmetric model such as the mass spectrum of superpartners 
and their interactions are highly dependent of how the SUSY is broken, which is also 
particularly of importance to search supersymmetric signals at collider experiments. 
Several breaking mechanism are known [3]. Among those, the gauge mediation scenario 
of the SUSY breaking |1] attracts attention, because it gives a solution to dangerous 
SUSY flavor problems which are generally caused by the introductions of superpart- 
ners of quarks and leptons. In the scenario, the breaking occurs at lower energy scale 
than those of other SUSY breaking models, so that the superpartner of the graviton, 
the gravitino, inevitably be the lightest supersymmetric particle (LSP). Here, we focus 
on the low-scale gauge mediation predicting the gravitino with O{10)eV mass in this 
letter. Such a scenario with very light gravitino is well motivated because it is com- 
pletely free from severe constraints coming from cosmology [5] such as the Big-Bang 
Nucleosynthesis [6] and the formation of large scale structure of our universe [7]. 

In the study of the gauge mediation scenario of the supersymmetry breaking, the 
gravitino mass ^3/2 is a very important parameter because m3/2 is directly related to 
the scale of the SUSY breaking, and also because the phenomenology and cosmology 
strongly depend on the gravitino mass. Thus, the experimental determination of the 
gravitino mass has a great impact on the study of the gauge mediation scenario. One 
way to determine the gravitino mass is to measure the decay width of superparticles 
into their superpartner and gravitino. In the gauge mediation scenario, the next lightest 
superparticle (NLSP) decays only into gravitino (and its superpartner), so the lifetime 
measurement of the NLSP gives a direct determination of the gravitino mass. 

Collider signals of the low-scale gauge mediation scenario depend on what the NLSP 
is. Though there are many candidates for the NLSP, we focus on the stau NLSP in this 
letter, which is predicted in wide parameter region of the scenario. Once produced, 
the stau NLSP eventually decays into a r-lepton and a gravitino with the lifetime of 
10~^^-10~^^sec if the gravitino mass is of O{10)eV, so that the decay length of the stau 
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NLSP is much shorter than the typical size of the detectors of coUider experiments. 
Thus, at the LHC experiment, for example, we expect the typical supersymmetric 
signal, namely, multi-jets associated with missing energy and r-leptons. Such a sig- 
nal is, however, generally predicted in various SUSY breaking scenarios. Thus, even 
though the existence of SUSY is likely to be confirmed by the LHC experiment, the 
measurement of the stau lifetime looks challenging at the LHC. 

The International Linear Collider (ILC) [8] is nothing but the experiment to pin 
down the new physics model under this kind of circumstance. The ideal environment 
with low background enables us to perform the precision study of the new physics. In 
this letter, we show that the measurement of the stau lifetime is possible by observing 
the distribution of the impact parameter, which is obtained by charged tracks caused 
by decay products of the r-lepton at the stau decay. In what follows, we consider how 
and how accurately the lifetime can be determined at the ILC based on a Monte-Carlo 
simulation study. 

2 Signal and Background Events 

The basic idea to measure the lifetime of the stau NLSP is the use of the simplest 
process producing the stau pair, e~^e~ — )■ f+f~, which eventually gives signal events 
involving two tau candidates with large impact parameters. In our analysis, we use 
only the events in which r-leptons decay hadronically in order to avoid the SM back- 
ground from the process e'^e^ — > W^W~ followed by the decay of the W boson into 
electron or muon (and neutrino). In addition, we assume that the contamination of 
the misidentified r-jet-like object into the r-jet sample is negligible. 

The cross section of the signal process, e+e~ f'^f", is depicted in the left panel 
of Fig. [Has a function of the stau mass nif. It can be seen that enough number of stau 
pairs are produced at the ILC. Stau NLSP decays into a r-lepton and a gravitino once 
it is produced, and its lifetime is determined by the formula, 

where Mpi ~ 2.44 x lO^^GeV is the reduced Planck mass, while m^i2 is the gravitino 
mass. When the gravitino mass is of the order of lOeV and the stau is lighter than 
200GeV, the decay length of the stau turns out to be longer than that of r-lepton 
[cTr — 87.11/im |9]) as shown in the right panel of Fig. [H 

Since the stau NLSP decays before reaching the first layer of the inner detector, 
which is now designed to locate at 16mm away from the beam pipe pi)], the charged 
track of the stau cannot be seen. Instead, we focus on the distribution of the impact 
parameter, which is the shortest distance to the track from the interaction point, of 
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Figure 1: Cross section of the signal process, e 
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the charged tracks arising from the stau decayi^^l Using the distribution of the impact 
parameter, we are able to measure the lifetime of the stau NLSP, which leads to the 
determination of the gravitino mass and therefore the scale of the SUSY breaking. 

If the decay length of the stau NLSP is much longer than that of the r-lepton, SM 
backgrounds are easily eliminated by requiring r-jets with large impact parameters. 
The SM backgrounds become, on the contrary, important when the decay length of the 
stau NLSP is the same order or shorter than that of the r-lepton. 

We expect large amount of SM background events including two r-jets. In the 
study of the impact-parameter distribution of the r-jets from the f decay, significant 
sources of SM backgrounds are expected to be coming from the processes, 

(a) [rr-BG] e+e" r+r" (+ ISR), 

(b) [lyiy-BG] e+e- ^ W+W- t+t-vv (+ ISR), 

(c) [ZZ-BG] e+e- ^ Z^Z^ t+t-uu (+ ISR), 

where ISR denotes the initial state radiation. Cross sections of these backgrounds are 
1144fb, 99fb, and 6fb for rr-BG, WW-BG, and ZZ-BG, respectively, when the center 



^^In our study, we neglect the effect of the magnetic field inside the detector, so that all the tracks of 
charged particles are approximated to be straight. Even if the charged tracks are not exactly straight 
in the actual experimental circumstance, the lifetime of the stau NLSP can be measured at the ILC as 
far as the profile of the magnetic field inside the detector is well understood. More complete analysis 
based on a full simulation study will be given elsewhere TT. 
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of mass energy at the collision is -y/i = SOOGeV. Notice that the cross sections for the 
WW- and ZZ-BG include relevant leptonic branching ratios of weak bosons. 

It should be noted that r-pair may be also produced by the two photon event 
77 — > r+r^, where photons are emitted from the initial-state e^. This type of 
background is, however, eliminated if an initial-state is detected; with the cur- 
rently expected performance of the ILC detector, the forward region is covered up to 
|?7| < 5.22 [I2], which will help to reduce the background. In addition, the two photon 
event is suppressed as we impose a lower bound on the energy of r+r~ system. In our 
analysis, we impose a lower bound on the total visible energy in the event, and neglect 
the background from the two photon event. Following Ref. [13], we require that -Evis, 
which is defined by the total energy of charged particles in the event, should be larger 
than 70GeV. (See the Cut 2 in the next section.) 

Now we discuss how we can eliminate rr-BG, VTiy-BG, and ZZ-BG. One important 
observation in reducing backgrounds is that, in the pair production processes (a) - (c), 
r^, W^, and Z° are likely to be produced with high rapidity. Then, in the lyiy-BG 
and ZZ-BG, are produced by the decay of weak bosons. Even in such cases, the 
information about the direction of weak bosons are inherited in the direction of decay 
products, i.e., r^. Thus, concentrating on events with r-clusters with low rapidity, 
standard-model backgrounds are significantly reduced. (See the Cut 3.) Because the 
cross section for the process e~^e~ — > r+r^ (+ ISR) is significantly larger than the 
signal cross section, the rapidity cut on r-clusters is not enough to reduce rr-BG. 
Large amount of remaining rr-BG can be eliminated by using the fact that, in the 
rr-BG, two r-leptons are (almost) back-to-back. If significant amount of momentum 
is carried away by ISR photons, and r~ are not back-to-back in general. However, 
even in such events, the ISR photon is likely to be (almost) parallel to the initial beam 
direction. Thus, we reject events in which two r-clusters projected on the x-y plane 
are back-to-back. (See the Cut 4.) We found that such cut is also useful to reduce the 
VrVT-BG. In addition, if an ISR photon is emitted with large transverse momentum, 
momenta of two r-leptons and that of the photon should be on one plane in the rr-BG 
event. This fact can be also used to reduce the rr-BG. (See the Cut 5.) 

3 Simulation Results 

In this section, we present several results of our simulation study. We have generated 
both the signal and background events using the HERWIG package [Mj. This package 
enables us to generate the events not only with momentum information but also those 
of decay points. In our study, we fix the stau mass to be nif = 120GeV. We neglect 
the uncertainty in the stau mass because the stau mass is expected to be determined 
with the accuracy less than 1% by measuring the energy distribution of reconstructed 
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Signal rr-BG WW-BG ZZ-BG 



Cut 1+2 

Cut 1+2+3 
Cut 1+2+3+4 
Cut 1+2+3+4+5 



2630 45172 3911 247 

2197 1798 911 79 

1316 525 362 26 

1307 116 353 26 



Table 1: The cut statistics for background suppression. Here, the number of events is 
normalized to the integrated luminosity C =100fb~^ with the center of mass energy y/s = 
500GeV. Mass of the stau NLSP is set to be 120GeV. In our analysis, the efficiency of the 
signal event is insensitive to Tf as far as the decay length of the stau is much shorter than 
the detector size, which is the case in the present setup. 

r-leptons from the decay of the stau NLSP [15] at the signal process. 

For each events, we identify decay products of r-leptons. The efficiency of the r-jet 
tagging is assumed to be 100%. For r-jet labeled by /, which is defined by hadronic 
decay products of I-th r-lepton, we define the momentum of the r-jet (Pj) to be 



where Pj is the three momentum of the particle i in the J-th r-jet, and the summation 
is over all the long-lived charged mesons in J-th r-jet. P/ is parameterized as 



Then, based on the discussion given in the previous section, we impose the following 
kinematical cuts on the generated events to reduce the SM backgrounds: 

1. The number of r-jet should be equal to 2. 

2. ^vis > 70GeV. 

3. I cos 01 1 < 0.85 and | costal < 0.85. 

4. -0.95 < cos(0i - 02) < 0.25. 

5. If there exist isolated photons with transverse momentum larger than 20GeV, we 
require |Pi x P2 x p^^''| > 0.1. (Here, p^^"* is the momentum of highest pt isolated 
photon, and a is the unit vector parallel to a.) 

In order to see how well the SM backgrounds are eliminated, and hence to see how 
well we can determine Tf, by imposing the kinematical cuts, we show the cut-table (the 
cut-statistics for background reduction) in Table [H It can be seen that the number of 
background events can be made small enough by imposing the kinematical cuts. 




(3.2) 



P/ = (P/ sin 9i cos (pi, Pj sin 9i sin 0/, Pj 



cos 9i). 



(3.3) 
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Figure 2: Distribution of the impact parameter for the signal and background events. The 
solid histograms are distributions of signals for Tf = 10~^^ (red) and 10~^^sec (blue), while 
the background is given by the shaded one. The stau mass is fixed to be 120 GeV. 



With the use of the generated events which pass the above kinematical cuts 1-5, 
we calculate the impact parameter of each r-jet, which is given by 



bi 



(3.4) 



where x/ is the decay point of the J-th r-lepton. Note that the summation over the 
index / should not be taken here. We expect that, in the realistic circumstance of the 
ILC experiment, the distribution of the impact parameter bj is obtained by measuring 
the shortest distance to the r-jet track from the interaction point. 

The resultant distributions of the impact parameter for the signal and background 
events are shown in Fig. |2j We adopt the error of the impact-parameter measurement 
of 5/xm here and hereafter according to Ref. [10]. It can be seen that the typical size 
of the impact parameter of the background events are ~ lOO/im, which is expected 
from the order of the decay length of the r-lepton. On the other hand, the typical 
size of the impact parameter for the signal event depends strongly on the lifetime f. 
This fact indicates that the lifetime of the stau NLSP is determined accurately once 
the distribution of the impact parameter is obtained accurately. 

Now, we quantitatively estimate how accurately the lifetime can be determined by 
using the distribution of the impact parameter. The number of event in i-th bin of the 
distribution is denoted as Ni^Tf); in our analysis, we have used the bins for < 6 < 
2mm with the bin-width of 20/im, so that the number of bins is 100. In order to see 
how well the lifetime can be constrained, we also calculate the theoretical expectation 
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Figure 3: The x^-variable of the £ = lOOfb ^ data for the three cases of the stau hfetime; 
3 X lO'^'^sec, 3 X 10~^^sec, and 3 x 10~^^sec. Stau mass is fixed to be 120 GeV. 



of the number of events in the i-th bin A^^ ( function of Tf. Here, A^^ ' {Tf) 

has been obtained by using the events generated with very high statistics. Then, for a 
given underlying value of the stau lifetime , we calculate the x^-variable as a function 
of the postulated value of the lifetime (denoted as rf"'^'^'^): 

4'""") - TssrrW E - A'f '(4'"""))' . (3^5) 

i Vf / i 

where i is summed over all the bins for the distribution of the impact-parameter. 

Taking C = lOOfb"^, in Fig. [31 we have shown the variable for three cases of 
the underlying stau lifetime; 3 x 10~^^sec (red line), 3 x lO^^^sec (green line), and 
3 X lO^^^sec (blue line). The minimum of the value divided by the degree of the 
freedom (100 — 1 = 99) is 0.93, 0.87, and 0.99, respectively. On the other hand, since 
Ax^{= — Xmin) = 1 corresponds to the observed lifetime with 68% C.L., the figure is 
showing that the lifetime of the stau NLSP can be measured to be (3.36tQ;g2) x 10~^'^sec, 
(3.13^°:°^) X 10-i3ggc^ (3.161°:°^) X 10-^2ggp^ respectively. 

We also consider the measurement accuracy of the lifetime of the stau NLSP. For 
this purpose, we estimate the Ax^ by replacing Ni{Tf) — )■ Nf^\Tr) in Eq. (13. 5p . Based 
on this quantity, we derive 68% and 95% confidence levels of upper and lower bounds 
on the lifetime. The result is shown in Fig. |H It can be seen that the nonzero lifetime 
of the stau NLSP can be confirmed when the lifetime is longer than ~ lO^^^sec. The 
measurement accuracy becomes better when the lifetime becomes longer. For instance, 
the lifetime can be measured with the accuracy less than 5% when > lO^^^sec. 
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Figure 4: Measurement accuracies of the lifetime of the stau NLSP. Upper and lower bound 
on the hfetime at 68% (red lines) and 95% (blue lines) confidence levels are shown as a 
function of the lifetime. The stau mass is fixed to be 120 GeV. 



4 Summary and discussions 

We have studied a collider signal with a stable gravitino of O{10)eV mass at the ILC. In 
our analysis, we particularly focus on the case that the NLSP is stau, which eventually 
decays into a gravitino and a r-lepton with the lifetime of 10^^^ — lO^^^sec. It has been 
found that the signal of the stau decay can be detected by measuring the distribution 
of the impact parameter, which is obtained by observing charged tracks caused by 
decay products of the r-lepton. We have also found that the lifetime can be measured 
with enough accuracy when it is longer than ~ 10~^^sec and the stau mass is about 
100 GeV. This fact guarantees the capability of the ILC to test the low-scale gauge 
mediation scenario of the SUSY breaking. 

In our study, we only consider the distribution of the impact parameter at the ILC 
with the center of mass energy 500 GeV to measure the lifetime of the stau NLSP. It 
may be possible to increase the significance to detect the signal by choosing appropriate 
center of mass energy according to the stau mass. In addition, the use of the incident 
polarized electron/positron beam will help us to reduce the background significantly. 
Furthermore, observing not only the impact parameter but also the direct information 
about the decay point of the r-lepton using 3-prong hadronically decay may enable us 
to measure the lifetime smaller than ~ 10~^^sec. These kinds of analysis require us to 
go to a full simulation of the detector performance and it is beyond the scope of this 
letter though those will be given in near future |11] . 
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Finally, we comment on the stau NLSP at the LHC. In the model of a long-lived 
stau, we also expect sizable number of SUSY events which contain r-jets with large 
impact parameters even at the LHC experiment. Although the full reconstruction of 
the SUSY event is challenging at the LHC, the determination of the lifetime of the stau 
NLSP may be possible. More interestingly, at the early stage of the LHC experiment, 
the discovery reach of the SUSY may be extended by looking for r-jets with large 
impact parameter. These subjects will be discussed elsewhere [16]. 
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